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Abstract The y-ray binaries LS 5039 and LS I +6T303 have been detected by 
Cerenkov telescopes at TeV energies, exhibiting periodic behavior correlated with 
the orbital period. These y-ray binary systems have also been recently detected by 
the Fermi Gamma-ray Telescope at GeV energies, and combination of GeV and 
TeV observations are providing both, expected and surprising results. We summa- 
rize these results, also considering the multi-frequency scenario, from the perspec- 
tive of pulsar systems. We discuss similarities and differences of models in which 
pulsar wind/star wind shocks, or pulsar wind zone processes lead to particles accel- 
erated enough to emit TeV photons. We discuss in detail the caveats of the current 
observations for detecting either accretion lines or pulsations from these objects. We 
also comment on the possibility for understanding the GeV to TeV emission from 
these binaries with a 2-components contribution to their spectrum. We show that it 
would be possible to accommodate both, normal pulsar emission and GeV / TeV 
fluxes that vary with orbital phase. We point out several aspects of this idea that are 
subject to test with data being currently taken. 



1 Prologue 

Four massive binaries have been discovered as variable very-high-energy (VHE) y- 
ray sources. They are PSR B1259-63 (Aharonian et al. 2005a), LS 5039 (Aharonian 
et al. 2005b, 2006), LS I +61°303 (Albert et al. 2006, 2009, Acciari 2008,2009), and 
Cyg X-l (Albert et al. 2007). HESS J0632+057 was found by the H.E.S.S. experi- 
ment (Aharonian et al. 2007) as one of the very few point like sources. Its positional 
association with the massive star MWC 148 led to the suggestion that its nature is 



Diego F. Torres 

Institucio Catalana de Recerca i Estudis Avancats (ICREA), Catalunya, Spain; and Institut de 
Ciencies de l'Espai (IEEC-CSIC), Campus UAB, Fac. de Ciencies, Torre C5, parell, 2a planta 
08193 Barcelona, Spain. e-mail; [dtorres@ieec . uab . es 1 



1 



2 



Diego F. Torres 



that of a gamma-ray binary system, albeit this has not yet been confirmed (Hinton et 
al. 2009). PSR B 1259-63 is obviously a system formed with a pulsar whereas Cyg 
X-l is most likely formed with a black hole compact object. The nature of the two 
remaining systems (LS I +61°303, and LS 5039, recently detected by Fermi, see 
below) is not yet settled. Variable gamma-ray emission was also reported from Cyg 
X-3, but it is only detected at GeV energies (Abdo et al. 2009c, Aleksic et al. 2010). 
The phenomenology presented by the systems that have been detected both at TeV 
and GeV energies can be distinguished in two classes: there are those presenting re- 
current TeV emission correlated with the orbit (the case of LS I +61°303, LS 5039, 
and PSR B1259-63) and there is also one case (Cyg X-l) that was hinted at in a flar- 
ing episode that has not been found to repeat yet. Cyg X-l and the three other TeV 
sources also differ in their SEDs. In Cyg X-l, the transient VHE luminosity was 
less than 1% of the X-ray luminosity. This Chapter analyzes some aspects of the 
high-energy and multi-wavelength phenomenology, especially of LS I +61°303 and 
LS 5039, in what concerns to their possible interpretation as non-accreting pulsar 
systems. 

In this same volume, contributions by R. Dubois et al. (Fermi collaboration) and 
J. Cortina give account of the observations at high and very-high energies of these 
sources. We shall only make notice of the spectral energy distribution, referring 
to the work of the former authors for details. Figure Q] shows the high and VHE 
spectrum for both LS 5039 and LS I +61°303. For LS I +61°303, the data plotted 
at the two energy bands are not corresponding to the same part of the orbit, with 
the Fermi data being an average. In the case of LS 5039, two broad band spectra 
-corresponding to the inferior and superior conjunction phase intervals- are shown 
for both energy bands, and although not contemporaneous, they cover multiple or- 
bital periods. 



2 Why is a non-accreting pulsar system a tenable alternative? 

In this section, we put some important aspects of the multi-wavelength information 
on LS 5039 and LS I +61°303 in context of the highest energy detection. We list the 
main reasons by which one can in principle sustain that a non-accreting pulsar sys- 
tem is a tenable alternative for these binaries (see also Dubus 2006, and Zdziarski, 
Neronov & Chernyakova 2008). 

1 ) The mass function does not provide definite constraints to the mass of the com- 
pact object, with a pulsar being an acceptable (if not preferred) alternative: 

Combining Newton's laws of gravitation and motion, the mass function is 
f(in\,m x ) = (47T 2 /G)((asin/) 3 /P^) = (m x sin;) 3 / (nt\ +m x ) 2 , where m\ and m x are 
the masses of the star and companion respectively, G is Newton's gravitational con- 
stant and i is the inclination angle of the binary orbit (defined so ; = 90° is edge 
on). The mass function has then units of mass, and is the minimum mass of the 
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Fig. 1 Summary of GeV to TeV observations of the two y-ray binaries detected with Fermi. High 
and very-high energy spectral data points of LS I +61°303 (left) and LS 5039 (right). Left: Fermi 
(integrating the whole orbit) data points are triangles; MAGIC data are lighter circles (high state 
phases 0.5-0.7); and VERITAS data are black circles (0.5-0.8). Right: Circles show the spectral 
data at INFC (TeV maximum, GeV minimum, phases 0.45-0.9; triangles show the spectral data at 
SUPC (TeV minimum, GeV maximum, phases < 0.45 and > 0.9). Higher energy data comes from 
H.E.S.S. Note that the data from the different telescopes are not contemporaneous, though they do 
cover multiple orbital periods. 



x-companion corresponding to a zero-mass m\ = and inclination angle ; = 90 : 
m x > f(mi,m x ). This means that the determination of the inclination and compan- 
ion mass is mandatory to determine the type of the compact object. If m\ = 10—15 
M©, and f(mi,m x ) ~ 0.01 the companion would be a neutron star if i > 25°. The 
mass function values found for LS 5039 and LS I +61°303 (see e.g., Casares et al. 
2005a,b; or Aragona et al. 2009) are low (~ 0.01). To constrain the compact ob- 
ject mass from the mass function, the orbital inclination and mass of the optical 
primary must be known, which are subject to their own uncertainties. The parame- 
ter space derived by Casares et al. (2005a,b), by Grundstrom et al. (2007), and by 
Aragona et al (2009) all allow for both a neutron star or a black hole in LS 5039 and 
LS I +61°303 the inclination of the orbit being poorly constrained (limits are not 
strict) 10° < i < 60, or even higher. It is interesting to note that the assumption of a 
neutron star in both objects, of 1 .4 M would imply a large inclination: the compact 
object would be a neutron star if i > 25° and a black hole otherwise. If ignorance 
of the orbital inclination can be represented by a random distribution of ;, the likeli- 
hood for the system to be formed by a neutron star is greater than if it is formed by 
a black hole compact object. 

2) There is no clear sign of an on-going accretion process: 

The X-ray and radio properties of LS 5039 and LS I +61°303 distinguish them 
from other X-ray binaries, and particularly for LS I +61°303, set it apart from other 
Be X-ray binaries (XRBs) (see below), since they both present no conclusive signs 
of accretion. If these systems accrete, given the sizes of their orbits, accretion would 
have to be wind-fed. Periodic ellipsoidal variations are also a common phenomenon 
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among high-mass X-ray binaries (HMXBs), provided that the optical companion 
fills its Roche lobe, due to tidal distorsions of the optical star. These are also not 
observed (Marti et al. 2004). The Bondi mass accretion rate is M ~ M w /(2r a /d s ) 2 
where M w is the stellar wind mass loss rate (typically about 10 7 ~ 8 M fl yr -1 ), d s is 
the orbital separation, and r a = 2GM c /v 2 v is the Bondi capture radius, with M c the 
mass of the compact object and v w the wind speed. This implies average values, 
e.g., for LS 5039 and assuming a polar flow of 2000 km s _1 , of 5 x 10 14 g s _1 . This 
value of mass accretion, when transformed into power, is close to, or even lower 
than the GeV luminosity of the system, implying unrealistic efficiencies (see, e.g., 
Dubus 2006). In order to avoid disc or magnetospheric accretion, the relativistic 
wind of the pulsar must be able to stop the infall of stellar matter, implying an upper 
limit to its putative period. Writing the spin-down power as a function of the pulsar 
magnetic field B and period P, it results P < 230B 1 / 2 M 15 '^ 4 ms, with M15 being the 
scale M/10 15 g s . Thus, if a fast millisecond pulsar is present in these systems, 
accretion onto it is not expected. 

As e.g., Zhang et al. (2010) confirmed with the analysis of deep INTEGRAL ob- 
servations of LS I +61°303, there is no high-energy cut-off at energies below 100 
keV, as it would be common in conventional accretion scenarios. Indeed, if the sys- 
tem is an accreting neutron star or black hole, one expects to find a cut-off power- 
law spectrum in the hard X-ray band with a cut-off energy normally at 10-60 keV 
for neutron stars (e.g., Filippova et al. 2005) and at ~T00 keV for black holes (Mc- 
Clintock & Remillard 2003). See also the LS 5039 Suzaku observations by Taka- 
hashi et al. 2009, where the X-ray spectral data up to 70 keV are described by a hard 
power-law with a phase-dependent photon index which varies within 1.45-1.61). 
Regarding spectral lines due to accretion, RXTE spectra of LS 5039 did appear to 
show a strong, broad Fe line (Ribo et al., 1999), although it was not confirmed by 
other X-ray missions (see, e.g., the XMM-Newton observations by Martocchia et al., 
2005). Several authors argued that being LS 5039 in the Galactic Plane, the Fe line 
is very likely due to Ridge emission sampled by the large RXTE field of view (e.g., 
Bosch-Ramon et al., 2005; Dubus 2006a, Zdziarski, Neronov, Chernyakova 2008). 
We discuss the prospects for line detection from these systems in greater detail be- 
low, where we put forward a note of caution given the relatively limited observation 
times that have been granted compared with the amount of time needed to detect 
accretion lines from these systems under some assumptions. 

3) If there is a pulsar, non-detection of radio pulses is expected: 

Zdziarski, Neronov & Chernyakova (2008) have already shown in detail why 
radio pulses would be absorbed in these compact binaries. We follow their derivation 
in what follows: The free-free absorption coefficient due to ions with the atomic 
charge, Z, is given by 



Off = (2V2^1/2 e 6 )/(33 /2 m 3/2 c)(/trr 3/2 z 2 WeWzV - 2 - ; 



(1) 
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where v is the frequency, «z is the Z-ion density, n e is the electron density, and g is 
the average Gaunt factorQ Using Z = 1, T = 10 5 K and v = 5 GHz, and averaging 
over Z, ~ 0.127 , - 3 / 2 (/i?/^ e )n?V _2 cm" 1 ~ O.nSr^nfv^crrr 1 , where T 
and v are in units of K and Hz, respectively, and jU e = 2/(l+X')~ 1.2 is the mean 
electron molecular weight. For the equatorial disc, the density can be expressed as 
n d,i{D) — «d,o(D / R-k)~^ \ where «d.o ~ 10 13 cm -3 , J — 3.2, and D is the distance 
from the center of the Be star (see, e.g., Waters et al. 1988, and the discussion be- 
low). The optical depth perpendicular to the disc plane is subject to a disc thickness 
of 2D tan Oq, with 6q the half-inclination angle. Thus, 
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where 3 x 10 12 cm corresponds to the periastron separation of LS I +61°303. Thus, 
the equatorial disc is optically thick to radio emission, including of course, pul- 
sations. This is consistent with the suppression of pulsed radio emission observed 
close to periastron in the system PSR B 1259-63 (Johnston et al. 1992, Melatos et al. 
1995). To calculate the radial optical depth of the fast polar wind from infinity down 
to a given value of D we can consider a clumpy wind (free-free absorption would 
still provide a high opacity even for smooth polar flows) where the wind density 
inside the clumps is 1 // times that of the smooth wind, («;) ~ n\f and 
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In a clumpy medium, % is an integral over the square of the density within the 
clumps times /, i.e., nff (or equivalently, an integral over (n;) 2 / -1 ). The opacity 
then results, 
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(4) 



Thus, according to the orbital solutions found for the systems, the neutron star would 
move in the optically thick region, explaining the absence of observed radio pulsa- 
tions. In addition of the high opacity, one can always entertain the possibility for the 
radio emission cone to be emitted off our line of sight, thus missing it entirely as 
in the case of any undetected-in-radio, bright GeV pulsars, even if at some portion 
of the orbit, radio pulses could escape. Pulsed emission in the X-ray or GeV band 
would be as important as those in radio to pinpoint a pulsar component, and we 



The Gaunt factor for hv <C kT and v v p (where v p is the plasma frequency) equals g = 
{3 1 / 2 /7t)[\n((2kT) 3 / 2 )/(7te 2 Zm e v)) - 5^/2], where 7e ~ 0.5772 is Euler's constant. 
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come back to these possibilities below. 

4) There are no clear signs of jets and the radio morphology is consistent with a 
pulsar component: 

From ~50 mas resolution radio images of LS I +61°303 obtained with Multi- 
Element Radio Linked Interferometer (MERLIN), extended, apparently precessing, 
radio emitting structures at angular extensions of 0.01 — 0.05 arcsec have been re- 
ported by Massi et al. (2001, 2004). However, recent Very Long Baseline Array 
(VLBA) imaging obtained by Dhawan et al. (2006) over a full orbit of LS I +61°303 
has shown the radio emission to come from angular scales smaller than about 7 mas 
(which is a projected size of 14 AU at an assumed distance of 2 kpc). This radio 
emission appeared cometary-like, and interpreted to be pointing away from the high 
mass star and thus being the smoking gun of a pulsar wind. A careful analysis of the 
images and phases, however shows that this interpretation is not straightforward: 
The tail is not always seemingly pointing in the right direction, at least when the 
polar flow of the Be star is considered. Nevertheless, there are some unambiguous 
results of these observations: a) No large features or high-velocity flows were noted 
in any of the observing days, which implies at least its non-permanent nature, b) 
The changes within 3 hours were found to be insignificant, so the velocity of the 
outflow can not be much over 0.05c. 

The MAGIC collaboration also conducted a radio campaign (in concurrence with 
TeV observations) to test these results using again MERLIN in the UK, the Euro- 
pean VLBI Network (EVN), and the VLBA in the USA (Albert et al. 2008). In 
this campaign, radio observations at different angular resolutions were conducted 
at the same time. It is then interesting to focus on the day where all these facili- 
ties observed at the same time (October 25-26, 2006). The results obtained by radio 
imaging at different angular scales show that the size of the radio emitting region 
of LS I +61°303 is constrained to be below ^6 mas (^12 projected AU), and the 
presence of persistent jets above this scale is therefore excluded. As in the case of 
Dhawan et al. (2006), these observations have shown a radio-emitting region extend- 
ing east-southeast from the brighter, unresolved emitting core. The outflow velocity 
implied by these observations is ~ 0.1c. The comparison between Dhawan et al. 
(2006) and Albert et al. (2008) images at the same orbital phase (but obtained 10 or- 
bital cycles apart) show a high degree of similarity on both its morphology and flux, 
which suggests periodicity and stability of the physical processes involved in the 
radio emission. These would be hard to attain if, for instance, the radio emission is 
the result of the random interaction of a steady flow with wind clumps. On the other 
hand, if the radio emission is produced by a milli-arcsecond scale jet, the required 
stability and periodic behavior of such a jet in order to produce the same radio map 
across many orbital cycles would be difficult to reconcile with the non-persistent 
nature of a larger scale (~100 mas) relativistic jet. Figure [2] shows that no jet can be 
seen at any scale. 

A similarly changing milli-arcsecond radio morphology was found for LS 5039 
(Ribo et al. 2008). In fact, early detection of elongated asymmetric emission in high- 
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Fig. 2 Radio images of LS I +61 303 obtained on 25 October 2006 with the VLBA (panel a) and 
on 26 October 2006 with MERLIN (panel b), EVN (panel c), and VLBA (panel d). In panel a, 
the 8 GHz VLBA image on 2 February 2006 (grey scale) is also overlaid, convolved with the 5 
GHz VLBA beam. This date corresponds to the same phase of LS I +61 303 (0 or b ~0.62). In all 
cases, the origin of coordinates is set at the VLBA peak of brightness on 25 October 2006, and the 
contours are drawn at (3,3 \/3,9,...) times the off-source rms. Figure reproduced from Albert et al. 
(2008), with permission of A AS. 



resolution radio images obtained with VLBA and EVN was interpreted as evidence 
of its microquasar nature, and suggested that the source was persistently producing 
jets with mildly relativistic ejections with a velocity of ~ 0.15c (Paredes et al. 2000, 
2002). Recent analysis by Ribo et al. (2008) show that two images obtained five 
days apart (recall the orbital period of about 4 days in this system, implying that 
the observations occurred in phases = 0.43 — 0.51 and 0.71 — 0.79) present a 
changing morphology. There is a core component with a constant flux density, and 
an elongated emission with a position angle (PA) that changes by 12° ± 3°, with the 
source being nearly symmetric in the first run and asymmetric in the second one. 
Analysis of the changes in the position angles, the inferred velocities of the outflows, 
and the steadiness of the radio flux, all make for an unpreferred interpretation in the 
framework of a microquasar model. As the authors conclude, they more naturally 
point to a young non-accreting pulsar scenario where the different morphologies 
detected at different orbital phases are due to the change of the relative positions 
between a pulsar and the companion star along the orbit. 

Dubus (2006a) has earlier explored this morphology change in pulsar wind/stellar 
wind scenarios, predicting a periodic change in the direction and shape of the ex- 
tended radio morphology as well as in the peak position of the radio core, depending 
on the orbital phase. The key concept is that whereas on a small scale, the shocked 
material flows away from the binary on a straight path, following the direction 
given by the vector difference of the stellar wind and orbital speeds v w — v or b; on 
a larger scale, the material shears with the orbital motion becoming important at 
~ d n ~ ocP or b/27t (about 1 AU for LS 5039)0 These effects combine to provide 
non-trivially shaped radio-maps, similar in aspect to those observed. The maps at 
0=0-0.25 are similar to the radio-morphology observed by Paredes et al. (2000), see 

2 This scale results from a comparison between the relevant flow timescale Tfl ow = d/v to reach 
a distance d from the shock and the orbital motion timescale T 01 t = d s /v ^ which can vary a lot 
for a highly eccentric orbit, and having taken v equal to its asymptotic value ac (a is the ratio of 
magnetic to kinetic energy, o" = (B 2 j 'Alt) j \T 'nm e c 2 )) . Thus, the shape and location of the cometary 
outflow can ultimately be used to constrain the magnetization. 
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Dubus (2006a). 

5) A young age of the system is consistent with observations: 

If the compact object is a young and energetic millisecond pulsar (necessary 
to avoid accretion onto it), the system should then be similarly young. Ribo et al. 
(2002) have put an upper limit to the age of LS 5039 by tracing back the proper mo- 
tion of the system to the plane. LS 5039 is a runaway system, with a total systemic 
velocity of about 150 km s and a component perpendicular to the galactic plane 
larger than 100 km s~ 1 . This is probably the result of an acceleration obtained during 
the supernova event that created the compact object in this binary system. Nitrogen 
enrichment in the atmosphere of the companion star also suggests a young age in 
LS 5039 (McSwain et al. 2004; rotational-mixing could also explain the enrichment, 
Casares et al. 2005b, see Dubus 2006). Note too that radio searches of supernova 
remnants have also failed for PSR B 1259-63, which has a measured spin-down age 
of 3 x 10 5 years (Johnston et al. 1992). 

6) Pulsar systems would be consistent with population predictions of Be XRBs: 

The Be XRBs, like LS I +61°303 are the most numerous class of XRBs known. 
At present, 64 Be XRBs are known in the Galaxy, and in 42 the compact object 
was confirmed to be a neutron star (NS) by the presence of the X-ray pulsations; 
in not a single one of them, a black hole was confirmed (Belczynski & Ziolkowski 
2009). Recent population synthesis made by the latter authors are consistent with 
0-2 black hole Be XRB in the Galaxy, which can perhaps be understood as a result 
of an evolutionary effect. The rotation of Be stars could be achieved during a period 
of Roche-lobe overflow mass transfer from its initially more massive companion, 
and if it looses most of the mass in the early pre-SN phase process, becoming a He 
star with mass of only a few M , when exploding as a supernova it can only leave 
a neutron star behind (see e.g., Tauris & van den Heuvel 2006). One has to bear in 
mind, however, that most of these Be XRBs are accreting systems, hosting pulsars 
with long periods (older spin-down systems) and only a few are expected to have 
a young age at any given moment. Spectral-wise, then, e.g., LS I +61°303 is very 
different from most of these objects above 10 keV, with the latter being cut off in 
hard X-rays (typically with ~20 keV bremsstrahlung-like spectra) with only very 
weak emission above ~400 keV as measured, e.g., by INTEGRAL. 



3 Caveats in the search for X-ray spectral lines 

The apparent absence of spectral lines in the X-ray spectra of Te V binaries has been 
often used as a proof of their non-accreting nature. However, a word of caution is 
useful in this respect. The X-ray continuum spectra of accretion dominated HMXBs 
are often described by a power law with photon index a ~ 1 — 2 (modified at higher 
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energies by an exponential cutoff between ^30-100keV). A spectrum of this form 
can be produced by inverse Compton scattering of soft X-rays by hot electrons in 
the accretion column near the compact object, and a part of this emission is scattered 
by the stellar wind of the massive companion. This results in a further non-thermal 
spectral component, but with a different absorption column depending on the orbital 
phase of the system. Furthermore, in some HMXBs a soft excess at ~ 0. 1 — 2 keV is 
detected, very common in systems hosting pulsars (Hickox et al 2004). In luminous 
systems this soft component can be explained by reprocessing of hard X-rays from 
the neutron star by optically thick, accreting material. For less luminous sources the 
soft excess is probably due to other processes, e.g. emission from photoionised or 
collisionally heated diffuse gas or thermal emission from the surface of the neutron 
star. On top of the continuum model, several spectral lines are usually present in 
these systems, neutral and ionized, such as Fe, Si, Mg, Al, N, Ca, mainly produced 
in the stellar wind or in the accretion disk (if any) illuminated by the strong X-ray 
emission of the compact object. One clear example of such a system, showing all the 
above mentioned spectral components is the HMXB 4U 1700-37 (composed by an 
O type star in a ^4 days orbit around a compact object of unknown nature: see e.g. 
Clark et al. 2002; van der Meer et al. 2005; see e.g. Figure[5]l. Note that in 4U1700- 
37 as well as many other binaries, the presence of spectral lines is highly dependent 
on the continuum spectrum and on the orbital phase of the system. 

So far a detailed high-resolution spectral analysis has been missing for the two 
debated TeV binaries: LS 5039 and LS I +61°303. The best data available for these 
kind of studies come from the XMM-Newton satellite, thanks to its large collecting 
area, spectral resolution, as well as the availability of its grating spectrometer. How- 
ever, at a given orbital phase, only very short observations have been taken (mainly 
aiming at monitoring the continuum spectral variability over the orbits). In the past 
XMM-Newton spectra of, e.g, LS 5039, an EW ~ 60 eV is the current la limit on 
the presence of Fe Ka (and only in a small part of the orbit), while in the available 
Suzaku observations, the limit on the detection of lines is 40 eV (see Takahashi et al. 
2009). Furthermore, with the current short pointings (at a given orbital phase) there 
are not enough counts to use anyhow the high-resolution spectral capabilities of the 
grating cameras, making impossible to answer any question about the presence of 
narrow lines. 

To understand what does this mean, one can consider the binary 4U 1700-37, a 
very similar HMXB to LS 5039 (i.e., similar companion star and orbital period, but 
no TeV emission, and clearly accreting) located at ~ 1.5 kpc rather than 3 kpc. Fig. 
0shows the XMM-Newton spectrum of 4U 1700-37 during an eclipse (van der Meer 
et al. 2005), and how would these accretion lines look, assuming all are present, 
in the available data for LS 5039. As a result of this simulation, one can see that 
none of the lines present in 4U 1700-37 would be detectable by the current data. In 
Figure[3]we also show the limits for the presence of lines derived from a 95 ks long 
Chandra recent observation of LS I +61°303 (Rea et al. 2010). 
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Fig. 3 Left panel: XMM spectrum of the 4-day binary system 4U 1700-37 at 0. 15 of phase (this is 
a reanalysis of observations made by van der Meer et al. 2005). Middle panel: simulation of what 
XMM-Newton would have seen of the 4U 1700-37's lines with the available short observations of 
LS 5039. Right panel: sensitivity plot for the presence of spectral lines in a ~ 100 ks long Chandra 
observation (the right two panels are adapted from Rea et al. 2010). 



4 Caveats in the search for pulsations 

The detection of pulsations is the only unambiguous tracer for a secure determina- 
tion of the pulsar nature of the compact objects hosted in LS 5039 and LS I +61°303. 
Deep searches for pulsations have been performed in the radio band at several fre- 
quencies, with the hope of detecting a fast spinning radio pulsar as in the case of 
the other TeV binary PSRB1259-63 (Johnston et al. 1999, 2005). However, no ra- 
dio pulsation have detected so far from any of these two sources. As it has been 
discussed earlier, this is anyway not surprising. In particular, note that at periastron, 
PSRB 1259-63 does not show radio pulsations, and that its periastron (given the 
large orbit, 3.4 year period) has about the same dimension of the major axis of the 
orbit of LS I +61°303 and it is way larger than that of LS 5039's. On the other hand, 
searches for pulsations in the X-ray band have many more chances of success than 
in the radio band. In fact, the X-ray pulsar beam is usually larger than the radio one, 
and the strong companion wind does not influence much the X-ray pulsed emission 
if present. However, what limits the X-ray pulsation search is the pulsed fraction 
sensitivity that current instruments can achieve. 

Until not long ago, archival observations which could give reliable upper lim- 
its on pulsations for fast spinning pulsars (P < 100 ms) hosted in LS5039 and 
LS I +61°303 were not very constraining, coming mainly from RXTE and XMM- 
Newton observations. In particular the high background of these instruments (es- 
pecially RXTE) limited the pulsed fraction sensitivity of these observations. For 
LS I +61°303 the deepest pulsed fraction limit was derived from a 41 ks XMM- 
Newton observation (Sidoli et al. 2006), which was < 28% (if not otherwise speci- 
fied, all pulsed fraction limits here are calculated assuming a sinusoidal profile and 
are reported at 90% confidence level in the 2-10 keV energy range), in the 12- 
200 ms period ranged Similarly, for LS 5039 the deepest limits for the presence of a 



3 Note that the RXTE monitoring observations performed in 1996 gave an upper limit of only 
32% in the 1-200 ms range, in fact, even though SXTE's timing resolution and collecting area 
was much larger than XMM-Newton, the much higher background contamination is a killer for 
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fast pulsar were derived from a 50 ks RXTE observation performed in 2003, giving 
an upper limit of Pf < 30%. Very recently, to increase the pulsed fraction sensitivity 
reducing the background contamination, long Chandra observations have been per- 
formed for both LS I +61°303 and LS 5039 (Rea et al. 2010). The results of these 
observations in the case of LS I +61°303 are summarized in very deep limits on 
the presence of X-ray pulsations, with an average limit of < 10% in the 6 ms - 10 s 
period range (Rea et al. 2010). Isolated rotational-powered pulsars emitting pulsed 
X-rays are all characterized by X-ray pulsed fraction much larger than the limits 
derived by Chandra for LS I +61°303. However, the presence of the strong stellar 
wind of the companion, and in particular the shock between the possible pulsar and 
the stellar wind, is a source of strong X-ray emission. This unpulsed X-ray emission 
coming from the shock can be responsible in diluting the pulsed X-ray emission 
proper to the pulsar, which then might easily end up to represent only < 10% of the 
total X-rays emitted by the systems. Of course, geometry and beaming can also be 
claimed for missing the pulsations at all. 

The possible detection of pulses could also come from observations of the Fermi 
telescope. Blind search pulsations of faint y-ray sources require very long obser- 
vation times, which implies the calculation of very large fast Fourier transforms 
(FFT). In addition, the significant frequency derivatives (/) typical of the y-ray pul- 
sars require that the FFT be repeated many times over a scan of /. Specifically, in 
order to keep the signal power within a single bin of the FFT, the frequency step 
should be Af = 1 /T oos , where T b s is the total observation time, and the steps size 
of / would have to be no larger than 1 /Tg b , which means an enormous number of 
/ trials for observation periods as long as ~1 year. Atwood et al. (2006) proposed 
the time-differencing technique, that reduces the number of / trials with a modest 
reduction in sensitivity. This method is based on the application of the FFT on the 
differences of the photon arrival times, rather than on the time series itself. In order 
to be efficient, only the time differences shorter than a predefined time window T w , 
that is significantly shorter than the whole observation, are considered in the FFT. 
In this way the number of required steps in / is reduced by a factor T w /T OPS . But 
unless the binary parameters are exactly known, they should also be scanned (like 
the /) to search for a solution. 

The correction of the photon arrival time series from the solar system reference 
frame to the pulsar reference frame is mostly affected by the Roemer delay in the bi- 
nary system. To evaluate it, five binary parameters are necessary: the orbital period 
(Pb), the epoch of the periastron (TO), the argument of the periapsis («), the pro- 
jected semi-major axis (Al), and the eccentricity (e). For LS 5039 andLS I +61°303, 
the first two parameters are know with relatively good precision, but the uncertain- 
ties on Al, e, and a> are larger than a few percent (e.g, see Casares et al. 2005; worse 
for LS I +61°303, see Aragona et al. 2009 and references therein). Thus, without 



detecting weak signals, resulting in a larger pulsed fraction limit. For that observations, Harrison 
et al. (2000) claims a limiting pulsed fraction of ~6%. However they considered the total count 
rate without correcting for the cosmic and instrumental background, which if corrected increases 
substantially the upper limit on the detectable pulsed fraction. 
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better knowledge of the binary system parameters, the computational difficulties in 
running a blind search make the detection of y-ray pulsations very challenging. 



5 Notes on the theoretical models based on pulsar systems 

High energy emission from pulsar binaries have been subject of study for a long 
time (e.g., Bignami et al. 1977, Maraschi and Treves 1981; Protheroe and Stanev 
1987, Arons and Tavani 1993, 1994; Tavani & Arons 1997; Bednarek 1997; Kirk 
et al. 1999, Ball and Kirk 2000; Dubus 2006a,b; Sierpowska-Bartosik 2008a,b and 
others). In most of these studies high-energy primaries are assumed to be accelerated 
at the shock formed by the collision of the (sometimes, putative) pulsar and the 
massive star winds. This generic feature allows them to be classified as wind- wind or 
inter-winds models. In other works the initial injection is assumed to come directly 
from the pulsar (the interacting particle population can be a result of equilibrium 
between this injected distribution and the losses to which it is subject, just as in 
the case of shock-provided electron primaries; and can also be the result of top- 
of-magnetosphere or inner-pulsar-wind shocks). In this case, and when opacities to 
y-ray production are high due to geometrical and physical characteristics, the pulsar 
wind zone (PWZ) could be the main origin of the high energy radiation, and we can 
classify this set as intra-wind or PWZ models. The closer the binary, the more these 
two scenarios are expected to produce similar results. Almost an order of magnitude 
difference in opacities can be seen between both systems discussed, with those in 
LS 5039 being higher, corresponding to the much shorter orbital periodicity (e.g., 
see Sierpowska-Bartosik & Torres 2009 for a comparison plot). Put otherwise, it 
would imply that a) intra-wind processes can not be completely neglected for LS 
5039, since the opacities for electrons in there is high and cascades can happen b) 
inter-wind models can not be completely neglected for LS I +6 1°303 (for these kind 
of more elongated systems, the probability for cascading in the pulsar wind region 
are lower and most of the energy carried by e in the PWZ can then be released in 
the shock region, where the local magnetic field traps the pairs to produce photons 
via synchrotron and IC processes). These models do not consider either possible 
contributions of a proton component (which can generate y-rays through pp or py 
processes (see Chenyakova et al. 2006). These extra components can, as discussed 
by these authors, also generate additional variability in the electron injection. These 
models do not contain a relativistic MHD approach of the collision of winds, but 
are rather based on simplifying assumptions of the hydrodynamical balance (see 
Romero et al. 2007 and Zdziarski et al. 2008 for a discussion on this issue). 

Through the study of intra-winds models, y-ray astronomy opens a window to 
study the electron distribution and magnetization properties of the pulsar winds in 
binaries, a virgin territory otherwise. Assuming intra-wind models, one can already 
prove that a mono-energetic distribution of leptons in the putative pulsar wind is 
ruled out by observations of LS I +61°303 and LS 5039. Assuming power-law dis- 
tributions in the winds of pulsars seems justified and possible, although there is a 
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not an priori expectation of the normalization and slope for them, and they so be- 
come free parameters of the models, which may in addition be subject to orbital 
variability, perhaps especially in close binaries. Although more complete models 
can certainly reduce/erase the need for such a change in the slope of the injection of 
leptons, as discussed, this is a caveat in the current versions. 

In any of these flavors of models, orbital periodicity is naturally encompassed. 
Short and random timescale variability can additionally be expected on top of this 
periodic behavior, for instance, as a result of granularity or clumping in the stel- 
lar wind, or the appearance of random shocks in the inner wind of the pulsar, which 
may modify local conditions off the average analysis discussed. The anti-correlation 
between GeV and TeV results is naturally encompassed in both of these flavors too 
(this is resulting just from the opacities evolution along the orbit and the system's 
geometry, this is a general output of binary models where the radiation is produced 
close to star even in the case of black holes). IC emission is enhanced (reduced) 
when the highly relativistic electrons seen by the observer encounter the seed pho- 
tons head-on (rear-on) i.e. at superior (inferior) conjunction. Inversely, VHE ab- 
sorption due to pair production will be maximum (minimum) at superior (inferior) 
conjunction. Gamma rays emitted in the vicinity of the compact object with energies 
above the ^30 GeV threshold inevitably pair produce with stellar photons (see e.g. 
Protheroe & Stanev 1987). On the contrary, emission in the Fermi range is largely 
unaffected by absorption but can be affected by cascading of higher-energy pho- 
tons. The phases of minimum and maximum flux in Fermi and H.E.S.S., as well as 
the anti-correlation, are then consistent with these generic expectations, suggesting 
Inverse Compton scattering is the dominant radiative process above 100 MeV with 
the additional effect of pair production (affecting photons above 30 GeV) further 
modulating the y-ray phenomenology (Bednarek 2006, 2007; Sierpowska-Bartosik 
& Torres 2008ab; Dubus et al. 2008). 

All in all, these models seem to work reasonably well for the VHE results, but 
both flavors fail in predicting the cutoff found by Fermi at GeV energies. Browsing 
SEDs from nominal models published smoothly connect both energy regions and 
predict levels of GeV-fluxes in excess of what is observed by Fermi (in the range 
~10-100 GeV). The Fermi cutoff is an unexpected feature in published models and 
it constitutes an observational surprise. 



5.1 A perspective on the GeV cutoffs ofLS I +61° 303 and LS 5039 

The passage of a compact object through a dense equatorial disk, such as the one 
commonly found in Be stars, would crush a putative pulsar wind nebula closer to 
the neutron star, increasing synchrotron losses and introducing a strong dependence 
with orbital phase of the electron energy distribution (e.g., Dubus 2006b). This may 
explain the appearance of the cutoff in the Fermi spectrum of LS I +61°303. How- 
ever, there is no such disk around the 06.5V star in LS 5039. Thus, the existence 
of the cutoff in both systems argues against explanations related with the properties 




Fig. 4 Interpreting the Fermi and TeV results of y-ray binaries with the 2-components idea. See 
text for details. 



of the companion, and it seems to require that the radiative process, or the primary 
injected population, are different in the low and high energy domains (Abdo et al. 
2009b). 

The GeV variability is also intriguing. If the GeV energy cutoff is interpreted in 
the outer gap model of pulsars, it is determined by the balance between acceleration 
and losses to curvature radiation, with the y-ray emission being pulsed. The absence 
of pulsations in Fermi data is not constraining, given the difficulty of detecting faint 
pulsars directly in gamma-rays, especially when in binary systems. But even when 
detected as DC radiation, GeV magnetospheric emission due to curvature radiation 
would be produced near the pulsar, and has no obvious reason to be modulated with 
the orbital motion. In addition, the GeV-TeV anti-correlation would also remain 
unexplained. Despite of this, the orbital-average spectral results for LS I +61°303 
and those found at the SUPC of LS 5039 do resemble the spectral properties of 
pulsars. The typical Fermi pulsar emission (of which dozens have been measured) 
has a hard power-law spectrum with a photon index in the range 1.0-2.0 and all of 
them present an exponential cutoff, the typical value of which is at about 2 GeV (see 
e.g. the Fermi pulsar catalog; Abdo et al. 2010). 

In what follows, we entertain the possibility that the GeV emission detected by 
Fermi is produced by two components. One would be the magnetospheric GeV 
emission coming from a (putative) pulsar in the system. As discussed, this emis- 
sion is expected to be steady along the orbit (i.e., unaffected by the orbital motion, 
since its origin is within the pulsar's light cylinder), and pulsed. The other compo- 
nent would come from regions located farther away from the pulsar, i.e., from either 
the inter-wind region or from the wind zone or from both: this radiation would then 
be unpulsed, and it is naturally expected to vary with the orbital phase. For clarity, 
we will refer to this un-pulsed component as the wind contribution. 

Fig. |4] graphically shows this idea for the richest GeV and TeV sets of data, 
corresponding to the observations of LS 5039 at INFC (TeV maximum, GeV mini- 
mum, phases 0.45-0.9; Fig.|4]left panel) and SUPC (TeV minimum, GeV maximum, 
phases < 0.45 and > 0.9; Fig.|4]right panel). Data comes from Fermi and H.E.S.S. 
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in the GeV and TeV energy ranges, respectively (Abdo et al. 2009b, Aharonian et 
al. 2006). The green curve in Fig. [4]represents a typical Fermi pulsar spectrum (a 
power-law slope of r = — 1.8 and an exponential cutoff at 2 GeV), and is here sup- 
posed to be steady along the orbit, hence reproduced in both INFC and SUPC panels 

of Fig. a 

There is nothing special in the shape of the assumed pulsar spectrum when com- 
pared with all pulsars found by Fermi. However, the flux level of the green curve is 
here constrained by observations: On one hand, it cannot be much larger, since oth- 
erwise (being unaffected by the orbital motion) the INFC LS 5039 spectrum would 
have been found at a larger flux level. On the other hand, one could have assumed 
a lower-flux pulsar contribution, but in that case the pulsar magnetosphere would 
become more and more irrelevant the less contributing it is at all phase ranges along 
the orbit, and a yet to be developed theoretical model would need to account for 
the whole range of observations, from GeV to TeV, with no magnetosphere con- 
tribution. Hence, the most interesting testing scenario is whether it is possible (and 
useful to understand the observed phenomenology) to have comparably contributing 
pulsar and wind components when the GeV emission is maximum (equivalently, a 
pulsar domination when the GeV emission is minimum). 

The light shaded area in Fig. |H between 10-100 GeV, emphasizes the region 
where no data is available at the moment. The other shadowed regions in Fig. [4] 
represent the minimum and maximum of the wind component allowed by the GeV 
data, given the assumption of the pulsar contribution, at INFC and SUPC. At INFC 
(left panel, red shadow) the assumed pulsar spectrum is consistent with the data 
itself, hence the minimum allowed value of the wind component is compatible with 
zero. In both plots, the dashed lines represent the wind contribution spectra assumed 
to derive the red and blue solid lines. These solid lines are the sum up of the pulsar 
(green line) and the wind (dashed lines) components. 

One can see from Figure [4] that the appearance of spectral cutoff and the GeV 
variability can be accommodated in the framework of this idea. In particular, if 
these two components are emitting GeV photons such that, e.g., at the SUPC both 
contribute similarly to the GeV flux (see Figure |4]right panel), then: 

• The anti-correlation of GeV- TeV fluxes is naturally maintained, since it is a 
generic feature embedded in inverse Compton models describing the TeV fluxes 
and the pulsar only sums up to it an orbitally-steady contribution. Essentially: 
from inverse Compton models, the GeV emission is enhanced (reduced) when 
the highly relativistic electrons seen by the observer encounter the seed photons 
head-on (rear-on), e.g, see Boettcher & Dermer (2005), Bednarek (2007). 

• Fermi sees pulsar-like spectra which varies with the orbital phase as the resulting 
effect of a changing dominance of the two contributions along the orbit. For in- 
stance, at 10 GeV in INFC the wind-related contribution, responsible of the TeV 
flux, may naturally take over and be larger than that of the pulsar, leading to the 
disappearance of the cutoff. The left panel of Fig.|4]shows that if one sums up a 
pulsed and a wind contribution, the cutoff may naturally disappear (or get very 
large) in INFC. In the case of SUPC (blue lines), instead, one is summing up a 
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pulsed contribution with something that also must increase (because of the pre- 
dicted GeV-TeV anti-correlation found in all inverse Compton models, and the 
constrain of the last Fermi upper limit) towards low energies and thus the cutoff 
is maintained. A possible increase of the energy (or total disappearance) of the 
GeV measured exponential cutoff (i.e., why there is no cutoff in the INFC of LS 
5039; and why we see such a high cutoff in the integrated-along-the-orbit spec- 
trum of LS I +61°303) can thus be accommodated within this idea. Fermi maybe 
summing up components that in its range of energies are changing dominance 
from the inner pulsar to the outer wind contributions. 

• It would be even harder to detect the pulsation in GeV data (since the pulsed 
emission is only a fraction of the y-ray flux observed). However, this concept 
also implies that a pulsation search in Fermi data will be more successful in those 
phase ranges corresponding to the Fermi minimum (INFC), where the pulsed 
fraction is necessarily larger than at SUPC, due to the hardness of the spectrum 
measured at the neighboring energy band. 

• Pulsar models based on wind processes (intra or inter-wind) must underproduce 
the detected Fermi emission, since they lack the contribution made directly by 
the pulsar magnetosphere. 

We caveat though that we presented here only a conceptual exploration of this 
possibility, hence the relative level between the contributions can be altered while 
maintaining the same overall idea. We also emphasize that the shape of the curves 
plotted in Fig. [4] for the pulsar and wind contribution is very general. The wind 
contribution responsible for generating the TeV radiation should be below the upper 
limit put by Fermi at 10 GeV; and, particularly at SUPC, it should contribute more at 
lower than at higher energies (if based on inverse Compton, e.g., Bottcher & Dermer 
2005, Dubus 2006, Sierpowska-Bartosik & Torres 2008). 

This 2-components idea can be tested by future observations. In particular, fur- 
ther observations by Fermi can confirm or rule out this concept by investigating the 
yellow energy range (as reported in Fig. increasing the observing time. Another 
possibility to test the 2 component model is with a very deep INTEGRAL exposure 
of LS 5039. The current INTEGRAL fluxes, derived from 3 Ms of IBIS/ISGRI data 
are (3.54±2.30) x 10~ n erg cmT 2 s _I for INFC, and a flux upper limit for the 
SUPC phase interval of 1.45 x 10~ n erg cm" 2 s" 1 (90% conf. level) in the 25- 
200 keV energy band (Hoffmann et al. 2009). The pulsar spectrum is not expected 
to be variable with the orbit, and the slope in the hard X-ray should be the extrap- 
olation of that in the MeV band (before the cut-off, of course). We verified that the 
extrapolation of the assumed pulsar contribution (green line in Fig.|4]i to low ener- 
gies was in agreement with these results, and especially with the INTEGRAL upper 
limit at SUPC. More IBIS/ISGRI data needs to be collected to possibly constrain the 
spectral slope of the putative pulsar in LS 5039. 

In several recent conferences where results of the VERITAS array have been 
presented (e.g. Aliu et al. 2010, this volume), it was noted that LS I +61°303 could 
now be in a low TeV-state at the usual phases where it was detected at such ener- 
gies; whereas, it seems to be shining normally (or even at a higher level) in GeV. 
LS I +61°303 is a more complex source than LS 5039, and Be-wind variations can 
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introduce variability in the orbital profiles. The latter has been already found in X- 
rays (Torres et al. 2010). With data now at hand, and lacking orbitally-resolved GeV 
and TeV follow up, we cannot say more as to the impact of the 2-components idea 
in the case of LS I +61°303. 

Finally, we remark that we are not ruling out that a comprehensive model of 
the wind shock or wind zone contribution that is producing the TeV radiation may 
-on its own- be enough to describe the GeV phenomenology in detail. We can 
only safely say that the current models did not appropriately predict it, despite the 
level of detail they already show. On the other hand, it can only be pertinent to 
consider that if these sources are pulsar-composed (see, e. g., Dubus 2006a), their 
magnetospheric contribution could be comparable to that found in all other pulsars 
detected by the Fermi LAT. If so, this idea presents a natural setup to understand the 
observed variability in the GeV spectra of these sources. 



6 Epilogue 

At the moment of writing, new VHE facilities are in the process of design. The 
CTA/AGIS class would put together an array of dozens of IACTs in order to enhance 
the sensitivity (aimed to be improved by one order of magnitude at 1 TeV when com- 
pared with H.E.S.S. /VERITAS) and energy acceptance (from a few tens of GeV to 
100 TeV). One of the aspects that such facilities could study is the formation of rel- 
ativistic outflows from highly magnetized, rotating objects. For instance, data on LS 
5039, analyzed in the context of PWZ models already rule out that mono-energetic 
electrons are responsible for most of the emission. What else can future data tell 
us? How would that impact on current models for dissipation in pulsar winds (see, 
e.g., Jaroschek et al. 2008 and references therein)? It is not implausible that close 
systems may trigger different phenomenology within the PWZ, ultimately affecting 
particle acceleration there. Models of particle energization and dissipation in pulsar 
winds are currently made for isolated objects, and we lack knowledge on whether 
the inclusion of such objects in close binaries will affect the wind behavior or even 
the magnetosphere in appreciable ways. Could we gain knowledge on this using 
future short-timescales y-ray observations? Also among the possibilities for future 
instruments, it is also worth noticing that of the determination of the duty cycles 
of high-energy phenomena by using continued observations of key objects (such as 
Cyg X-l) with current instrument's sensitivity using sub-arrays of future ACTs. It is 
interesting to note that neutrino detection or non-detection with ICECUBE will also 
shed light on the nature of the y-ray emission / limit the neutrino-to-photon ratio 
irrespective of the system composition (e.g., see Aharonian et al. 2006b, and Torres 
& Halzen 2007 for applications to LS 5039 and LS I +61°303, respectively). For 
instance, Neronov & Ribordy (2009) presented a hadronic model for y-ray binaries 
in which the multi-TeV neutrino flux from the source can be much higher and/or 
harder than the detected TeV flux, and where most neutrinos are produced in pp 
interactions close to the bright massive star, in a region optically thick for the TeV 
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photons. The secondary pairs of these processes would participate in the emission 
at lower energies. The nature of LS 5039 and LS I +61°303 y-ray binary systems 
is still unknown. Whereas current data are consistent, and in some cases prefer, an 
interpretation based on non-accreting pulsars a final proof in either way (pulsations, 
accretion lines) is pending. 

This work has been supported by grants AYA2009 -07 391 and SGR2009-81 1. The 
author warmly acknowledges A. Caliandro, R. Dubois, G. Dubus, D. Hadasch, N. 
Rea, and A. Sierpowska-Bartosik for discussions. 
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